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Abstract: If the waterfall field of hybrid inflation couples to a U{1) gauge field, 
the waterfall can generate a statistically anisotropic contribution to the curvature 
perturbation. We investigate this possibility, generalising in several directions the 
seminal work of Yokoyama and Soda. The statistical anisotropy of the bispectrum 
could be detectable by PLANCK even if the statistical anisotropy of the spectrum 
is too small to detect. 
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1. Introduction 

Although there is so far no evidence for statistical anisotropy of the primordial cur- 
vature perturbation mechanisms have been proposed for generating it. Most of 
them invoke a vector field. 

One mechanism takes the vector field to be homogeneous during inflation, but 
causes significant anisotropy in the expansion |l[] (for a recent review of this approach 
see [0). Then the perturbations of scalar fields generated from the vacuum fluctu- 
ation will be statistically anisotropic, and so too will be ( on the usual assumption 
that it originates from one or more of these perturbations. 

We here invoke a different mechanism ^ (for the most recent paper on this 
approach see Taking the inflationary expansion to be practically isotropic, this 

*^Thc use of a vector field to generate a contribution to C was first mooted in [pj. 
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mechanism generates a perturbation of the vector field from the vacuum fiuctuation, 
which in turn generates a contribution to (. Results using this mechanism are at 
best approximate, because the unperturbed part of the vector field will cause at some 
level anisotropic expansion which generates additional statistical anisotropy through 
the first mechanism, but existing calculations ignore that effect and we will do the 
same. 

We work with the setup of The vector field is a U{1) gauge field coupled to 
the waterfall field of hybrid infiation.*^ The dominant contribution to ( is supposed 
to come from the perturbation of the infiaton field. But the perturbation of the 
gauge field is supposed to generate an additional contribution during the waterfall 
that ends infiation. The waterfall is taken to be practically instantaneous. We extend 
the original treatment of the scenario in several respects. First, we do not assume 
the infiaton potential V = Vo + m^0^/2 (which is ruled out by observation). Second, 
we do not assume that the perturbation of the gauge field is exactly scale-invariant. 
Third, we take into account the time-dependence of the gauge field. 

We take for granted the main ideas of modern cosmology described for instance 
in and use the notation and definitions of 0, §]• The unperturbed universe has 
the line element 

ds'^ = -dt^ + a^{t)5ijdx'dxK (1.1) 

For any smooth rotationally invariant quantity g{:s.,t), uniquely defined during some 
era, we can choose a slicing (fixed t) and threading (fixed x) and then write 

g{^,t) = g{t) + 6g{^,t). (1.2) 

Going to a different slicing with a time displacement 5t(x, t), we have to first order 

[6g{^, i)] - 6g{^, t) = g{^, t) - (?(x, t) ^ -g{t)bt{^, t). (1.3) 

We will invoke this 'gauge transformation' without comment. In most cases g is 
homogeneous on one of the slicings. 

2. The curvature perturbation ( 

2.1 Definition and 5N formula 

To define ( one smoothes the metric on a super-horizon scale, and adopts the co- 
moving threading and the slicing of uniform energy density p. Then |^ 

C(x,t) = (5[lna(x,t)] = (5[ln (a(x, t)/a(t))] = 5iV(x,t), (2.1) 

where a(x, t) is the locally defined scale factor (such that a comoving volume element 
is proportional to a^(x, t)) and a{t) is its unperturbed value. The number of e-folds 

*^The study of this setup with non-Abelian gauge fields can be found in [R. 
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of expansion A^(x, t, t^) starts from a slice at time on which a is unperturbed ('flat 
slice') and ends on a uniform p slice at time t. Since the expansion between two flat 
slices is uniform, 6N is independent of t^,. The change in ( between two times is 

C(x, t2) - C(x, ti) = (5iV(x, h, ts), (2.2) 

where now both the initial and final slices have uniform p. 

By virtue of the smoothing, the energy conservation equation is valid locally: 

pW = 3^^(p(t)+P(x,t)). (2.3) 

In consequence, C = during an era when P{p) is a unique function. The success 
of the BBN calculation shows that P = p/3 to high accuracy when cosmological 
scales start to enter the horizon. Then ( has a time-independent value ({x) that is 
strongly constrained by observation. Within observational errors it is gaussian and 
statistically isotropic. Its spectrum is nearly independent of k, with 

V(;{k) ~ (5 X 10-^)2. (2.4) 

For the reduced bispectrum |TD| /nl, current observation give |/nl| ^ 100 and barring 
a detection PLANCK will give |/nl| ^ 10. For /nl to ever be observable we need 
I/nlI > 1. 

We will work to first order in so that 

C(x,t) = i7(t)5V (2.5) 

where 6tfp is the time displacement from the fiat slice to the the uniform-p slice. 
A second-order calculation of ( is needed only to treat very small non-gaussianity 
corresponding to |/nl| ^ 1- 

To explain the near scale- invariance of the observed (fc), it is usually supposed 
that N{-x,t) is determined by the values of one or more fields 0i(x, t), evaluated 
during infiation at an epoch when relevant scales have left the horizon: 

Ar(x,t)=Ar(0*(x),0;(x),---,t). (2.6) 

The fields are defined on a fiat slice and denoting their values by 0* we write 

0:(x) = 0* + 50*(x) (2.7) 

and [ill 

C(x, t) = J2 NMK^) + N^Mi^)S<l^*M) + ■ ■ ■ ' (2-8) 

where a subscript i denotes d/(f)* evaluated at 0*(x) = 0*. The 0, are usually taken 
to be scalar fields, but it has been proposed 0, Q that some or all of them may be 
components of a vector field. 
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On each scale /c, the field perturbations are generated from the vacuum fiuctua- 
tion at horizon exit. Ignoring scales leaving the horizon after Eq. (|2.8| ) defines a 
classical quantity (, which is independent of the choice of t*. In general it depends 
on t, settling down to the observed quantity by some time tf. 

Cosmological scales have a fairly narrow range AA; ~ 15 or so. Choosing t^, as 
the epoch when the shortest scale leaves the horizon, scalar fields with the canonical 
kinetic term, that satisfy the slow-roll approximation, have a nearly Gaussian uncor- 
related perturbations with spectrum Vs<j)i* — (-ff/27r)^. To have the observed nearly 
gaussian C(x) Eq. ( |2.8|) has to be dominated by one or more linear terms (at least 
when t = t{). Keeping only linear terms, 

Pc(x,t):^^iV2n^. + ---, (2.9) 

where the terms exhibited correspond to scalar fields, and the dots indicate vector 
field contributions. The contribution of the latter is positive like the rest 0]. The 
non-linear terms may give non-gaussianity that can be observed in the future. 

2.2 Slow-roll inflation 

Infiation corresponds to e// < 1 where en = —H/H"^. During infiation each coordi- 
nate wavenumber k (scale) leaves the horizon when k = aH, and we are interested 
only in the era of infiation after horizon exit for the biggest observable scale k ~ ao-ffo- 
Here H = a(t)/a(t), k is the coordinate wavenumber and the subscript denotes 
the present. We need ej^ <C 1 at least while cosmological scales leave the horizon to 
generate the nearly scale-invariant V^{k). 

We are interested in single-field slow-roll infiation. Here, the only field with sig- 
nificant variation during infiation is the infiaton. Its unperturbed value 0(t) satisfies 
the slow-roll approximation. 

3H^ ~ -V'{(j)), (2.10) 
e ^ ^M^{V'/Vf ^ « 1 (2.11) 
p{t) = 3Ml,H^ ~ 1/(0). (2.12) 

The perturbation 50^, can be removed by a shift 5t^,(x, t), which means that it 
generates a time-independent contribution to Since the /nl generated by is 
negligible [|l^, we can work to first order in 50*, 

C0(x) = -(if/0)50,(x). (2.13) 

Choosing as the epoch of horizon exit, we find the spectral index is given by 

n^{k)-l = dVi;Jd\nk = 2ri-6e, (2.15) 
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where rj = MpV"/V with |?7| <ti 1, and the right hand sides are evaluated at the 
epoch of horizon exit aH = k. 

Although it is not our central concern, we mention at this point the case of 
multi-field slow-roll inflation, where two or more fields vary significantly. Taking 
to be the field pointing along the trajectory at horizon exit, Eq. ( p.l4| ) still applies 
to that case. 

If ( depends only on the part of the action that we are considering, C(p{k) can be 
identified with the observed quantity C(^)- More generally we have 

^^(5x10-^. (2.16) 

This inequality is important for two reasons. First, it makes the tensor fraction 
r < 16e. Then the slow roll approximation gives what has been called the Lyth 
bound, on the variation A0 of the infiaton field after the observable universe leaves 
the horizon. Without any assumption about the function e(0) after the first few 



e-folds, one finds 13 



10-1 (A0 /Mpf > 16e > r. (2.17) 

If e(0(t)) increases with time, 10^^ in the above expression is replaced fl^ by 0.0003. 
According to these results, an observable r cannot be obtained with A0 ^ Mp 
(small-field model). 

The other use of the inequality is for curvaton-type models, where ( is generated 
almost entirely by the perturbation of some field that has a negligible effect during 
inflation. Then will be negligible compared with V(, and r will be negligible 
compared with 16e so that it is unobservable. 

All of this assumes slow-roll inflation, in which it is assumed that there is no 
time-dependent field except the slowly-rolling infiaton fields. If there is such a field 
the shift in the initial time generated by 50* will be accompanied by a shift in the 
value of that field, which could allow to be time-dependent and avoid the inequality 



The possibility of avoiding this inequality was mooted in ||T5|, |T6l but they did not 



find a mechanism. One can easily avoid the inequality by abandoning the canonical 



kinetic term for the infiaton [|T7| and we are for the first time pointing to a possible 



mechanism with the canonical kinetic term. 

3. The model 

3.1 Hybrid inflation 

The relevant part of the action is 



S = d^x\f^^ 



\mIR - \d,<^d^<^ - \d,xd'x - \f\cl>)F,,F^^ - V 



(3.1) 
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with F^j, = dfj^Bij — dyB^ and 5^ a f/(l) gauge field. Following [Q we use the gauge 
with Bq = diBi = 0, and work with Ai = Bi/a which is the field defined with respect 
to the locally orthonormal basis (as opposed to Bi which is defined with respect to 
the coordinate basis). The raised component is = A^ (as opposed to -B* = Bi/a?). 
We also define W = /A, which would be the canonically normalized field if / were 
constant. To fix the normalization of /, we set / = 1 at a time just before the 
waterfall begins. 

Before the waterfall the potential is 

V{<P, X, A) = V, + V{cp) + im2(0, A)x^ + ^Ax' + \li^A^ (3.2) 
m'{(P,A) = h^A^ + g^"" -m". (3.3) 

The waterfall field x is supposed to be the radial part of a complex field which is 



charged under the t/(l), generating the first term of Eq. (|3.3| ). 

This is the usual hybrid infiation potential [|18], ^ except for the presence of 
A. We assume that the values of the parameters and fields give what has been 



called standard hybrid infiation |^0[. At each location, the waterfall begins when 
m^(0(x, t), A(x, if:)) falls to zero. Before it begins, the waterfall field x vanishes up 
to a vacuum fiuctuation which is set to zero, and we have slow-roll infiation with 

V = Vo + V{ct>) ~ V,. (3.4) 

Cosmological scales are supposed to leave the horizon before the waterfall begins. 
We will take H to be constant which is typically a good approximation. 

3.2 Field equations 

During the waterfall, x moves to it's vev and then infiation ends. We will assume 
that the duration of the waterfall is so short that it can be taken to occur on a 
practically unique slice of spacetime. The evolution of cj) and A is therefore required 
only before the waterfall begins. 

To work out the field equations, previous authors have taken /(0) to be a function 
of time with / oc a"(t), and have taken spacetime to be unperturbed. Then the action 
( p.l|) gives for the unperturbed fields 

0(t) + 3if0(t) + V'{<P{t)) = (3.5) 

W(t) + 3HW{t) + i2^W{t) = 0, (3.6) 

where 

jj^ = H\2 + a){l-a). (3.7) 

By virtue of the fiatness conditions on the potential (e ^ 1 and \ri\ ^ 1), the 
first expression is expected to give the slow-roll approximation ( p. 10 ) more or less 
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independently of the initial condition. Similarly, the second equation is expected to 
give the slow-roll approximation 

3HW ~ -/i^w (3.8) 

if l/ip ^ which we assume. 

In terms of W, the coupling h'^A'^x'^ becomes h'^W'^x^^ where h = hf. Before 
horizon exit on cosmological scales we are taking W to be a practically free field 
corresponding to h 1. With a ~ 1 this would give at t = t„ a tiny coupling 
h -C e~^'= which would have practically no effect. We therefore assume a ~ —2. 

For the first order perturbations, f (x a" gives 

6Mt) + SH6Mt) + {{k/af + V"(0(t))(50k = (3.9) 
'5Wk(t) + 3H6Wk{t) + {{k/af + /i^)5Wk = 0. (3.10) 

Since and W are slowly varying, this flat spacetime calculation is expected to 
hold in the perturbed universe on the flat slicing. It is expected because ||, ^] the 
effect of the metric perturbation (back-reaction) on Eq. (|3.10|) is proportional to the 
small quantity W(t). 

Since we are assuming /(0), the choice / oc a" corresponds to 

/ oc exp j^y2e{(t))M^]-^d(^ . (3.11) 
This gives the perturbation 

Since / is a function of 0, the term oc fF^^F^^ in the action couples cj) and W 
so that the right hand sides of Eqs. ( |2.10|) , ( p. 81) , (|3.9|) and ( p. 10]) are nonzero. We 
calculate them in the Appendix, and show that they are negligible if 

where pw is the energy density of W. We will assume this condition. Note that it 
implies pw ^ P, which is anyway needed because we are taking the expansion of the 
universe to be isotropic. From Eq. (|2.16|) , the condition is guaranteed if W/H < 10^. 



3.3 Spectrum of W 

The evolution equation for W(x, t) is the same as that of a free scalar field with 
mass-squared p?, and we are assuming <^ H^. Treating 5W\^ as an operator 
and assuming the vacuum state well before horizon exit gives the approximately 
scale-independent vacuum expectation value 

^(«l'^(O^H-i,(0> = - *-A^) ^»(k + k') (-) (-^) (3.14) 

n^cc - 1 = 2p^/3H^, (3.15) 
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where hats denote unit vectors. According to Eq. (3^), 5Wk has constant phase 
which means that it can be treated as a classical quantity with this correlator. 
We are interested in t = t„, when 



a{t„)H 



Also, we are interested in cosmological scales, which have a range AA^^ ~ 15 and a 
typical central value A^jt ~ 50. 
The decomposition 

W(x,t) = W(t) + 5W(x,t) (3.17) 

is made in some box of coordinate size L around the observable universe, with W(t) 
the average within the box. After smoothing on a cosmological scale k, the spatial 
average of {6WY (evaluated within a region not many orders of magnitude bigger 
than the observable universe) is of order ln{kL)H and we assume that the box is not 
too big, ln(/cL) roughly of order 1. Guided by these results, we assume W{t) ^ H, 
which is reasonable because W'^(t) at a typical position is expected to be at least of 
order the mean square of {6Wy evaluated within a much larger box p2|. 



4. Including the waterfall contribution 

4.1 End-of-infiation formula 

At an epoch t+ just after inflation ends, 

C(x,t+) =C+(X) =C<^(X)+Cw(x) 



(4.1) 



where Cw is the waterfall contribution. 

To calculate Cw, we suppose that the waterfall happens very quickly so that it 
can be regarded as taking place on a single spacetime slice. Then [BOl, |2r 



H5tuU) = H 



^Pw(x) _ 6p^ 

P{tw) 



H 



f){U 



H6t 



(4.2) 



In this equation, 5ti2(x) is the proper time elapsing between a uniform-p slice at time 
just before the waterfall, and a uniform-p slice at time t+ just after the waterfall. 
Because |p| is much smaller during inflation than afterwards, Sti2 is practically the 
same 6tp„, the displacement from the initial uniform-p slice to the waterfall slice. 
Using Eq. (p75|), we see that 

C(x,t+)=ff5t(x), (4.3) 

where 6t is the displacement from the fiat slice at t„ to the waterfall slice. 

This end-of-infiation formula actually holds if the waterfall slice is replaced by 
any sufficiently brief transition from inflation to non- inflation. In it is invoked 



- 8 - 



for the transition beginning during the waterfall, at the epoch when the evolution of 
X becomes non-linear. We are here applying it to the entire waterfall. It was first 
given |]23| with A in Eq. ( p.3| ) replaced by a scalar field. In the slope of the 



potential in the A direction was assumed to be negligible corresponding to single- 



field hybrid inflation, and the same assumption was made in several later papers [E4 



The assumption was relaxed in [E^, EBI, 12^, corresponding to what has been called 



multi-brid inflation. Following we are here taking A to be the magnitude of 
a f/(l) gauge field. One can also replace A by a non-Abelian gauge field 0, |28|, p9| . 

Before continuing, let us ask what is required to make the waterfall sufficiently 
brief. Just to have the error |A^„| ^ one presumably needs At(x) ^ |(5ti2(x)|. 



which is equivalent to 0].*'^ 

HAt < rll"^ < ^ 5 X 10-5. (4.4) 

During the waterfall m^{(j),A) goes from to — of the waterfall, with > H^. 
We therefore expect At to be at least of order 1/m, and we need m/H -C \/ M-pjH 
so that A ^ 1 Hence Eq. ( [4.4| ) requires an inflation scale if/Mp < 10^ GeV 

or 1/1/^ < IQi^GeV. This rules out GUT hybrid inflation iy^l^ ~ lO^^ or so) 
but easily allows inflation at the scale of supersymmetry breaking (if of order the 
gravitino mass < 10^ GeV or so). 

A stronger requirement might be needed to justify a calculation of the non- 
gaussianity parameter /nl, because it refers to the non-gaussian part of C that is only 
of order 'P^^'^^/nl ^ 10~^. A reasonable estimate for |A(^„/(^„| might be HAt/vl^^. 
Then, in the worst case that A(^w is completely non-gaussian, one would require the 
very low inflation scale H/Mp <C 10~^^/^l. We proceed on the assumption that the 
inflation scale is sufficiently small. 

4.2 Waterfall contribution in our model 

Without at first specifying the nature of A, we now calculate St. The fields on the 
waterfall slice have values given by m^(0„(x), A„(x)) = 0. It was noted in that 
the time dependence of a waterfall slice could be important. Thus let us define a 
'time- dependent waterfall slice' t) by 

m2(0„(x,t),A(x,t)) = 0. (4.5) 

(If this equation has more than one solution 0w(x, t) we choose one of them.) If 0„ is 
negligible, the waterfall occurs when 0(x, t) falls to the practically time-independent 
waterfall slice 0w(x). If instead (p is negligible, the waterfall occurs when the time- 
dependent waterfall slice 0„(x, t) meets the practically time-independent field value 



^■^The second inequality allows for contributions to C(x) that might be generated after inflation 
by fields different from (jj and A, and assumes that the contribution of the latter undergoes no 
further change. 
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0(x). Since we deal with hybrid inflation is negative, but (f)^ might have either 
sign. If it is also negative we need |0w| < \<P\ or inflation will never end. 

If (5t(x) is the displacement from the flat slice at t„ to the waterfall slice, this 
gives to first order va. 5t 

0(X, U + 5t(x)) = 0(tw) + 50(X, tw) + 0(tw)5t(x) 
0„(X, t„ + 5t(x)) = 0w(tw) + 50w(x, tw) + 0w(iw)5i(x), 

where 5<^w is defined on the flat slice. Setting 5t = gives 0(tw) = 0w(iv 
5t we have 

/-( + \ _ ZJ^^r \ _ ^^<Pw{^:U) -(50(x,tw) 



(4.6) 
(4.7) 

Evaluating 



i75t(x) = H- 



Now we invoke Eq. 



(x,t) = -(m^ 



C(x,t+) = C0(x) 1 + 



0(t„) - 0w(iw) 

). Discounting the strong cancellation 

Using Eqs. ( glSD , (pj), ( ^Tl and (p]) we get#^ 

1P 1 + 2X^2 J ^ " 



1 + 2XA^ ■ 



where C?i(x) is defined by Eq. (|2.13|) and 

Cw = -X (W(t^) ■ 5W(x,t„) + {5W{^,K)f) 
X = h^/y/2eMpmg. 



(4.8) 
/i^A^ it gives 
(4.9) 



(4.10) 

(4.11) 
(4.12) 



Previous authors except ignored the time-dependence of A, which means that 
they implicitely set 0w(iw) = to obtain = Cw*^ 

Taking = t„, the second term of Eq. ( [4.10D is the contribution of to (, 
which means that the first term is the contribution of 50*. It differs slightly from the 
result found earlier in Eq. ( [2.1ij| ), but the difference is not significant; to calculate the 
contribution of 50* taking account of W, one would have to include the anisotropy 
of the expansion of the universe caused by W which is beyond the scope of our 
investigation. Therefore, at the level of our calculation there is no change in the 
usual assumption that ^^(x, t) is constant. 

The formalism that we have given involves S(j)^, which is a function of SA and 
hence of both 6(f) and 6W. A more direct approach is to use 4>vi(t) defined by 
m^{(f)^, (iy//(0„))2). Then 50„ is a function only of W, and 

leading directly to Eq. ( |4.10|) . 

^'^Terms involving a product of 5(j) with itself or 5W are dropped because they are negligible. 

"^^^The spectrum of the waterfall contribution found in |^ is negligible (smaller than ours by a 
factor (1 + 2XA'^)'^e~'^^''). We have not been able to follow this calculation, which is not from first 
principles because A is treated as a scalar. 
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4.3 Anisotropic spectrum and bispectrum 

Since W ^ H, the linear term of Eq. ( [4.1 1]) dominates, leading to 



C 



1 -/3 ( A-kV 



(4.14) 



/3 = ^^-^^e^- (4-16) 



^ l + 2Xif2' ^^-'^^ 

where 

Current observation requires /3 < 10~^, and barring a detection PLANCK will give 
(3 < 10~^ Using Eqs. ( ^.4p and ( |2.14| ), the observed value of requires 

Including the second term of Eq. ( 4.11 ) we find ^ 

/NL = /^T(l + r"'(ki,k2,k3)) (4.18) 

where 

/- = HA-kO'^-(A-b)^+ (k. .fa)(A-k.)(A.fa) ^ ^ ^^^^^ 

E /^3 

(with ki + k2 + ka = 0) and 

^ f- - l±M^/3^ ~ 3 (100/3)^/^ ^ 1 + 2XA^ 

(For the final expression we used Eqs. ( p.l4|) , ([4.16| ), and ( [4.17|) .) The last two 
factors, omitted in the original calculation [^, allow PLANCK to detect /^l even if 
it does not detect (3. 

This calculation ignores the time-dependence of e. Allowing time-dependence for 
e would multiply Cw by a factor [e(tcos)/e(i^w)]^^^, where tcos is the epoch of horizon 
exit for a typical cosmological scale. The factor might be significantly different from 
1, but there is little point in including it because its effect is indistinguishable from 
the effect of the tilt factor e^. 



5. Conclusion 

If the waterfall of hybrid inflation is sufficiently brief, it takes place on a practically 
unique slice of spacetime. Then the waterfall slice contributes to the curvature 
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perturbation if its location depends on some field A different from both the infiaton 
and the waterfall field. 

We generalised in several directions the model of Yokoyama and Soda [Q, that 
takes A to he a U (1) gauge field. The model makes ( statistically anisotropic, and we 
find that the prediction for /nl could be verified by PLANCK, even if the prediction 
for the anisotropy of "P^ is too small to be detected. 

The weak point of the model is the special form Eq. that is required to 

get the gauge kinetic function /(0) oc a~^. We are not aware of any well-motivated 
hybrid infiaton potential that would lead to a well-motivated This is in contrast 

with the case of non-hybrid inflation [^, where one can take /(0) and V{(l)) to have 
exponentially increasing behaviour that might be reasonable in string theory [Q, 
and which could correspond to an attractor (late-time limit) pl| . 

In the course of our investigation we noticed that the presence, during slow-roll 
inflation, of a time-dependent field different from the infiaton might allow a significant 
decrease in the spectrum of the curvature perturbation after horizon exit. That does 
not however happen in our case. 
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A. Equations of Motion for 0(x, and W(x, 

Extremizing the action in Eq. with respect to fields 0, -B^ and their derivatives 
we obtain field equations 

[d, + d. In v^] + V'{<P) + ^ff'{<f>)F^,F^^ = 0; (A.l) 

[d, + d,\n^]fF^'' = 0, (A.2) 

where g = det^g^^^) and /^^ = df /d(f). Choosing the temporal gauge Bq = and a 
line element of the unperturbed universe in Eq. one finds equations of motion 

for the fields 0(x, t) and B(x, t) 

+ 3/^0 _ + V'{<P) = (A.3) 

Bi + (h + 2^] Bi- a'^V^Bi = a-^2^djBi, (A.4) 
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Recasting the above equations in terms of W = fB/a and dropping gradient terms, 
one arrives at equations of motion for homogeneous fields and W(t) 



f 



(A.5) 



W + 3HW + \ 2H^ - Hj - J \ W = Q, 



(A.6) 



where we also used ^ 0. 

Decomposing the field W(x, t) as in Eq. (|3.17|) and similarly the field (^(x, t), 
we find equations of motion for perturbations 50(x, t) and 5W(x, t) from Eqs. (|A.4| ) 
and ( |A.4| ). Keeping only the first order terms and switching to the Fourier space 
they become 



+ 3H6(f) + {k/af 6(j) + V"6(f) 



7 



5W + \ H - j] 5W - 6 Ij \ W 



,(A.7) 



5W + 3H5W + \ 2H^ - Hj - J \ 5W - a'^VHW 



w. 



(A.8) 



For exponentially varying gauge kinetic function / in Eq. the above expressions 

become 



+ 3H6(f) + {k/af 6(f) + V"6(j) 
2a 



/2eMp 
5W + 3H6W + i2^5W 
aW 



W + H{l-a)W 

2 



6W + H{l-a)6W + a^W6(j) 



,(A.9) 



-1 SW 

a 



2eMp 



6(1) + H {l + 2a)6(f)+ {-] 6(f) 



(A.IO) 



where /x^ = (2 + a){l — a)H'^ and a ^ —2. 

The energy density of the vector field in Eq. (|3.1|) is given by pb(x, t) 



—pF^yF^^ /A. From this it is easy to see that the background value of pBi^,t) is 
given by 



W+ \H-^\W 



(A.ll) 
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The right hand side of Eq.( A.5 ) is neghgible if ps satisfies Eq. ( |3.13 ). We now show 
that the same is true of the right hand sides of Eqs. ( |A.9|) and ([A.10|) . At the epoch 
k ~ aH, the terms on the left hand sides are of order and Eq. (|3.13|) ensures that 



the right hand sides are indeed much smaller. At the epoch aH/k = exp{Nk(t)) ^ 1, 
the first term of each left hand side is negligible. The other two terms are of order 
\r]\ = \ V"\/3H^ for Eq. and of order \r]w\ = \p^\/5H^ for Eq. ( [ATTUD . Eq. (^) 

ensures that the right hand side of Eq. ( |A.9| ) is negligible, and it ensures that the right 
hand side of Eq. (|A.10|) is negligible if also \riw\ ^ 10^^. But the latter condition is 
irrelevant, because its violation makes the time-dependence of W (coming then from 
the right hand side) negligible. 
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